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ABSTRACT 

We present a reanalysis of the relationship between asteroid albedo and polarization 
properties using the albedos derived from the Wide-field Infrared Survey Explorer. We 
find that the function that best describes this relation is a three-dimensional linear fit 
in the space of log (albedo)-log (polarization slope)-log (minimum polarization). When 
projected to two dimensions the parameters of the fit are consistent with those found in 
previous work. We also define p* as the quantity of maximal polarization variation when 
compared with albedo and present the best fitting albedo-p* relation. Some asteroid 
taxonomic types stand out in this three-dimensional space, notably the E, B, and M 
Tholen types, while others cluster in clumps coincident with the S- and C-complex 
bodies. We note that both low albedo and small (D < 30 km) asteroids are under- 
represented in the polarimetric sample, and we encourage future polarimetric surveys 
to focus on these bodies. 

1. Introduction 

As light scatters off the surface of atmosphereless bodies, it is instilled with a small linear 
polarization. The degree of linear polarization of the scattered light measured by the observer is 
a function of the phase angle of observation and the composition and structure of the surface, in 
particular the interrelated parameters of albedo, index of refraction, and space between scattering 
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elements (e.g. lMuinonenlll989l ; IShkuratov et a/.lll994l ). Early work quantified the relation between 
phase angle (the angle between the direct ion to the sun and the observer as seen from the target, a) 
and polarization ( Dollfus &: Zellnerll979 ) and this effect can be used in parallel with the magnitude - 
phase effect to probe the scattering physics of atmosphereless surfaces (jMuinonen et aZ.ll2002l . 120091 ) . 



As expected from classical scattering models, the light reflected from a surface is polarized 
perpendicular to the scattering plane for large phase angles, which is referred to as a positive 
polarization. For small phase angles, however, light acquires a polarization in the scattering plane 
due to an increase in the dominance of second-order scattering. This case is referred to as negative 
polarization, as it is perpendicular to the positive case and thus carries a negative sign when the 
polarization coordinate system is rotated to account for the viewing geometry. The angle where 
the phase curve transitions from positive to negative is referred to as the inversion angle («o)- By 
definition the value of the polarization must go to zero at a = 0°, though some work has suggested 
that s urfaces may have a sec o ndary trough at very small angles related to the optical opposition 
effect (IRosenbush et aLlll997l ). ICellino et all (|2005al ) find no evidence for a polarimetric opposition 
effect in their sample, though high albedo objects are not represented there. 

From the studies of the scattering properties of the lunar surface, a relationship was found 
between albedo and the p arameters used to describe the polari metric-phase effect of the lunar 



regolith (jBowell et aZ.lll973l ) that was then extended to asteroids (jZellner et aZ.lll974j ; ICellino et al. 
19991 ). of the form: 



\ogpv = C\ log h + C 2 
log p v = C 3 log P min + C 4 



(1) 

(2) 



where py is the geometric albedo, h is the linear slope of the phase curve at the inversion angle, 
and P m in is the value of the largest negative polarization (i.e. the depth of the negative trough), 
usually expressed as an absolute value. We show an illustration of these parameters and two typical 
polarization-phase curves in Figure [H We note that the polarization shown in this figure is the P r 
value that has been rotated to account for viewing geometry, such that P r > is the amplitude 
of polarization perpendicular to the scattering plane and P r < is the amplitude parallel to the 
scattering plane. A polarization component ±45° from the scattering plane is typically not observed 
at any phase angle for asteroids, and thus is ignored in this diagram. 



Cellino et all (|1999l ) present the most recent best-fitting values for the constants in the above 
equations: d = -1.118±0.071, C 2 = -1.779±0.062, C 3 = -1.357±0.140, and C 4 = -0.858±0.030. 
In this work, we revise the best-fitting values for t hese constants in lig ht of new albedo data from 
the Wide-field Infrare d Survey Explorer ( WISE, IWright et al\ |2010| ) and the planetary science 
extension NEOWISE (|Mainzer et aZ.ll2011al ) . The aim of this work is two-fold: firstly, while WISE 
provides us with albedos for a large fraction of the known asteroids, calibration of this relationship 
will allow it to be applied to objects that were not observed by WISE; secondly, the behavior of the 
polarization of asteroid surfaces helps us determine the surface mineralogy and this relationship 
represents a critical component in this determination. Through application of both thermal infrared 
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and polarimetric data we can gain a broader understanding of the behavior of asteroids across the 
Solar system. 




Fig. 1. — Illustration of the typical polarization-phase behavior for two different types of asteroids, 
with the inversion angle (an), the minimum polarization (P m i n ), and the slope of polarization at 
the inversion angle (h) labeled for the dashed blue curve. Example curves for a generic S-type 
asteroid (solid red line) and a generic C-type (dashed blue line) are shown. 



2. Data 



We draw our list of polarimetric properties for asteroids from a range of sources. The domi- 
nant con tributor is the Astronomica l Polarimetric Database presented in the Planetary Data System 
(PDS) (ILupishko &: Vasilyevl 120081 ) which was a compilation of the polarimetric properties of in- 
dividual asteroids in the literature up to t he data of pub l ication. We als o incorporate va l ues fo r 
h, Pmin, and/or a n for asteroids given by: Celling et al\ (11999I. l2005al a): iFornasier et al\ (12 006); 



20091 ) : lBelskava et all (|2010|) 



Gil-Huttonl (|2007ah : lGU-Hutton et all (|2007bl . l2008h : lMasiero k Cellinol 
We note that as these data are drawn from a range of different instruments, uncertainties in the 
absolute calibration may result in a larger scatter than is actually present. A comprehensive survey 
of polarimetric properties of a large number of objects conducted with a single instrument would 
reduce this possible source of error, and so is strongly encouraged. 

Determination of h, P m i n , and «o au require polarimetric measurements spanning a range of 
phase angles. The inversion angle can typically be determined to a reasonable level of accuracy 
with a few bounding measurements at a ~ 20°. Polarimetric slope is more difficult to determine, 
especially for objects located farther from the sun that are rarely observable at phase angles much 
beyond the inversion angle (e.g. objects that do not come within ~ 2.9 AU of the sun can never 
be observed at phase angles a > 20°). Careful timing of observations can ensure adequate phase 
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coverage that will allow for an accurate determination of the slope. The depth of minimum polar- 
ization is often the most difficult parameter to determine for some objects, as it requires observing 
at small phase angles that are not frequently available for asteroids in the inner Main Belt, Hun- 
garia, Mars Crosser, and NEO populations. Additionally, determining this value requires evenly 
spaced observations over the full branch of negative polarization, rather than just a few bounding 
measurements as required for both «o an d h. As such, relative errors on P m %n tend to be larger 
than measured for the other polarimetric parameters. Where errors on polarimetric parameters 
were not given by the source, we assume values based on the errors on the published data in those 
sources. 

The albedos we use for this work are drawn from the values derived for Main Belt asteroids 



(MB As) published in lMasiero et al\ (120111 ). F or objects in the NEO or Mars Crossing populations, 
we draw albedos from the appropriate lists ( Mainzer et al. 2011e f|) which were derived using a 
method identical to that used for the MBAs. All of the objects with both defined polarimetric 
phase curves as well as WISE-determined albedos had low identifying numbers, implying that they 
were some of the first objects discovered, and thus likely to preferentially sample the largest minor 
bodies of the solar system. These large asteroids were more likely t o have been se e n in m ultiple 
bands by WISE which allows for fitting of the beaming parameter. iMainzer et al\ (|2011bl ) show 
that in cases such as this the error on albedo as an absolute measurement is ~ 20% of the measured 
albedo value, however internal comparisons are better than this limit. 

A primary concern in any analysis of albedos derived from infrared-determined diameters is 
the quality of the optical measurements used. We draw our H magnitudes from the Minor Planet 
Center's orbital element catalog (MPCORE0), as discussed in iMasiero et al\ (|201ll ) . While other 
studies have found an offset between measured magnitudes and thos e predicted from the H absolute 



magn it ude value, with objec ts being fainter than expected (e.g. lJuric et al. 



2002 



Parker et al 



20081 ) . IMainzer et al\ (|2011bl ) find that, in general, no offset corrections to H are required for 
the most recent releases of the magnit udes. An exception t o this result has been found for some 



(120111 ) and IMainzer et al\ (j201ld ): many of 



objects with unusually high albedos in lMasiero et al. 
these object s are coincident in orbital-element-space with the Hungarias and the Vesta family. 
Harris et al. (jl989l ) found that the commonly assumed value of G = 0.15 is inappropriate for these 
types of high-albedo objects, and a value of G ~ 0.4 may be more appropriate. Revising G to 
this value would result in an offset of up to ~ 0.3 mag in the H magnitude depending on the 
initial H fit, however this correction is not required for most objects. In the past, photometric 
measurements for many asteroids that contributed to the H magnitudes in the MPCORB catalog 
were acquired with unfiltered CCDs. New, filtered observations and refined handling of previous 
photometry have largely mitigated the effect of unfiltered measurements on H values. (T. Spahr, 
2012, private communication). 



Mainzer et al\ (j201ld ) present a comparison of the WISE albedos to the IRAS albedos and 



http://www. rninorplanetcenter. net/iau/MPCORB. html 
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find a good match for most objects, though some scatter is seen, especially at the smallest sizes 
where the IRAS signal-to-noise ratios were poor compared to WISE. This albedo error assumes 
moderate-to-low light curve amplitudes and well characterized H and G values. This error will 
result in an uncertainty in the offsets of the linear fits (i.e. C2 and C4 in Equations [1] and [2]) though 
the slopes should be u naffected. We include this error in our fits below. We note that recently 
Muinonen et all pOlCj) have introduced a three-parameter photometric system (H,Gi,G2) to better 
characterize the behavior of the photometric phase effect which may reduce some of these errors, 
but we note that this system requires accurate photometry over a large phase window, which is not 
available for many asteroids. 

We show the polarimetric and albedo data used for this work in Table Q] (ellipses indicate an 
unmeasured polarimetric property). When these two data sets are combined, we have 65 objects 
with measured albedos and polarimetric slopes, and 112 with albedos and minimum pola r izatio n 
values. This result is a improvement compared to the data set presented bv lCellino et al\ (|1999l ). 
who performed a similar analysis using IRAS albedos of 37 objects for the slope-albedo fit and 
16 for the minimum polarization-albedo fit. We note that due to the brightness requirements of 
most polarimeters and the polarimetric survey strategies employed, these lists are dominated by the 
largest known asteroids. Approximately half our sample have sizes over 100 km, and three-quarters 
are larger than 50 km. Thus while the largest asteroids are well sampled, there is a distinct lack of 
small bodies in these li sts. We also note th at despite the fact that low-albedo objects dominate the 
Main Belt population (jMasiero et adl201ll ). they are under-represented in the polarimetric surveys 



(see below). As WISE is sensitive to thermal infrared light, the detection p robability for asteroid s 
is effectively unbiased with respect to the albedos of the objects observed (jMainzer et a/.ll2011el ). 
and thus the distribution of albedos seen with WISE is a more accurate representation of the true 
population than is the distribution seen for optically-selected samples. Extending polarimetric 
coverage to both smaller sizes and low albedo objects through a large-scale campaign is critical to 
extending and generalizing the trends seen here and in previous work. 



Table 1:: Compiled asteroid albedos and polarimetric prop- 
erties 



Asteroid 


Pv 


a 


h 


p 

1 mm 


2 


0.142 ±0.018 


18.1 ±0.1 


0.228 ± 0.003 


1.38 ±0.05 


5 


0.245 ±0.051 


19.1 ±0.1 


0.096 ± 0.050 


0.70 ± 0.05 


6 


0.269 ± 0.049 


20.8 ±0.2 


0.091 ± 0.050 


0.80 ± 0.05 


8 


0.261 ± 0.048 


20.0 ±0.1 


0.104 ±0.003 


0.68 ±0.05 


9 


0.134 ±0.016 


21.8 ±0.1 


0.102 ±0.003 


0.74 ± 0.05 


10 


0.058 ± 0.005 






1.50 ±0.05 


11 


0.158 ±0.036 


18.9 ±0.2 


0.124 ±0.003 


0.73 ± 0.05 


12 


0.140 ±0.014 


20.8 ±0.2 


0.121 ±0.003 


0.73 ±0.01 


13 


0.069 ± 0.022 


21.7 ±0.5 


0.257 ±0.003 


2.10 ±0.05 
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Table 1:: (continued) 



Asteroid 


Pv 


a 


h 


1 mm 


14 


0.221 ±0.022 


20.5 ±0.2 


0.105 ±0.003 


0.82 ±0.10 


15 


0.206 ± 0.055 


20.6 ±0.2 


0.087 ±0.005 


0.72 ± 0.02 


17 


0.160 ±0.009 




0.131 ±0.003 


0.74 ± 0.05 


18 


0.221 ± 0.082 


21.6 ±0.1 


0.101 ±0.003 


0.87 ±0.05 


19 


0.050 ± 0.020 


21.7 ±0.2 


0.305 ± 0.003 


1.72 ±0.05 


22 


0.169 ±0.061 






0.83 ± 0.04 


24 


0.064 ± 0.016 




0.191 ±0.003 


1.63 ±0.10 


27 


0.201 ± 0.058 




0.099 ± 0.003 


0.70 ± 0.05 


29 


0.157 ±0.035 


22.0 ±0.2 


0.098 ± 0.003 


0.88 ±0.10 


30 


0.171 ±0.034 


19.8 ±0.5 


0.104 ±0.003 


0.78 ± 0.05 


31 


0.045 ± 0.044 






1.32 ±0.10 


32 


0.230 ± 0.065 






0.63 ±0.05 


39 


0.245 ± 0.056 


21.0 ±0.1 


0.090 ± 0.003 


0.79 ± 0.05 


40 


0.195 ±0.019 


20.8 ±0.2 


0.100 ±0.003 


0.85 ±0.05 


46 


0.052 ±0.011 






1.54 ±0.10 


47 


0.067 ± 0.009 




0.204 ±0.003 


1.44 ±0.05 


51 


0.100 ±0.026 




0.292 ±0.050 


1.86 ±0.05 


54 


0.049 ± 0.008 


22.2 ±0.5 


0.357 ±0.050 


1.95 ±0.05 


56 


0.050 ± 0.006 


19.7 ±0.2 


0.318 ±0.003 


1.47 ±0.05 


57 


0.182 ±0.047 






0.71 ±0.10 


58 


0.059 ± 0.005 






1.70 ±0.10 


63 


0.159 ±0.028 


19.8 ±0.1 


0.102 ±0.003 


0.70 ± 0.05 


64 


0.676 ± 0.223 


18.2 ±0.2 


0.036 ±0.001 


0.32 ±0.05 


68 


0.207 ± 0.025 






0.68 ±0.05 


70 


0.040 ± 0.009 






1.83 ±0.05 


71 


0.248 ± 0.035 




0.061 ±0.005 


0.61 ±0.02 


73 


0.186 ±0.018 






0.76 ± 0.06 


75 


0.098 ± 0.014 


20.2 ±0.1 


0.103 ±0.017 




77 


0.153 ±0.046 






1.25 ±0.14 


80 


0.182 ±0.026 






0.75 ± 0.05 


83 


0.086 ± 0.021 






1.47 ±0.10 


84 


0.053 ±0.017 


20.3 ±0.5 


0.306 ± 0.050 


1.49 ±0.05 


85 


0.063 ± 0.025 






1.36 ±0.10 


89 


0.185 ±0.034 




0.119 ±0.050 


0.90 ±0.05 


95 


0.056 ± 0.009 






1.78 ±0.05 


97 


0.206 ± 0.046 


22.1 ±0.1 


0.174 ±0.018 




113 


0.223 ± 0.031 




0.081 ±0.005 
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Table 1:: (continued) 



Asteroid 


PV 


a 


h 


p . 

1 mm 


114 


0.088 ± 0.010 






1.24 ±0.10 


115 


0.654 ±0.124 






0.71 ±0.05 


118 


0.139 ±0.031 






0.80 ±0.12 


121 


0.077 ±0.010 






1.72 ±0.05 


125 


0.115 ±0.027 




0.145 ±0.033 


0.83 ±0.02 


129 


0.157 ±0.026 






0.90 ±0.05 


131 


0.164 ±0.011 




0.208 ±0.059 




132 


0.120 ±0.008 




0.146 ±0.006 


1.13 ±0.03 


135 


0.152 ±0.050 






1.06 ±0.10 


138 


0.161 ±0.028 




0.103 ±0.020 




139 


0.045 ± 0.023 




0.262 ±0.050 


1.31 ±0.05 


141 


0.049 ± 0.010 


20.6 ±0.5 


0.330 ± 0.050 


1.78 ±0.05 


145 


0.043 ± 0.004 






1.86 ±0.05 


153 


0.046 ± 0.008 






1.05 ±0.05 


182 


0.210 ±0.059 






0.64 ±0.05 


184 


0.107 ±0.019 






0.93 ±0.06 


188 


0.157 ±0.055 




0.140 ±0.015 




189 


0.199 ±0.024 






1.26 ±0.10 


192 


0.288 ± 0.040 


19.8 ±0.1 


0.084 ± 0.003 


0.75 ±0.05 


197 


0.239 ± 0.026 






0.79 ±0.08 


201 


0.097 ± 0.006 






1.00 ±0.05 


204 


0.163 ±0.044 






0.83 ±0.12 


216 


0.111 ±0.034 






1.27 ±0.05 


217 


0.044 ± 0.005 






0.82 ±0.05 


230 


0.171 ±0.076 


20.6 ±0.2 


0.122 ±0.003 


0.94 ±0.05 


234 


0.151 ±0.034 


27.0 ±2.0 




1.60 ±0.20 


250 


0.112 ±0.021 






0.88 ±0.08 


259 


0.042 ± 0.005 






1.25 ±0.05 


270 


0.254 ± 0.043 






0.65 ±0.05 


305 


0.182 ±0.028 






0.64 ±0.10 


306 


0.174 ±0.060 






0.66 ±0.10 


324 


0.063 ± 0.012 


20.0 ±0.1 


0.278 ± 0.003 


1.46 ±0.05 


334 


0.051 ±0.016 






1.32 ±0.05 


338 


0.165 ±0.028 






0.98 ±0.10 


345 


0.059 ± 0.012 






1.55 ±0.05 


347 


0.213 ±0.041 


22.6 ±0.1 


0.113 ±0.003 


0.78 ±0.03 


349 


0.153 ±0.018 






0.39 ±0.05 



-8- 



Table 1:: (continued) 



Asteroid 


PV 


a 


h 


1 mm 


351 


0.171 ±0.046 






0.74 ± 0.09 


354 


0.173 ±0.032 






0.51 ±0.10 


356 


0.053 ± 0.015 






1.50 ±0.10 


377 


0.056 ± 0.025 


19.8 ±0.2 


0.206 ±0.005 


1.76 ±0.04 


384 


0.190 ±0.040 






0.94 ±0.35 


396 


0.139 ±0.025 






1.34 ±0.09 


409 


0.050 ±0.010 


19.9 ±0.2 


0.191 ±0.005 




410 


0.043 ± 0.007 




0.313 ±0.050 


1.94 ±0.05 


415 


0.086 ± 0.009 






1.28 ±0.10 


423 


0.066 ± 0.005 






1.40 ±0.05 


441 


0.139 ±0.026 






1.41 ±0.12 


451 


0.069 ± 0.006 






1.62 ±0.05 


466 


0.086 ± 0.009 






1.60 ±0.10 


511 


0.073 ± 0.006 


19.4 ±0.1 


0.277 ±0.003 


1.69 ±0.05 


532 


0.202 ± 0.039 




0.122 ±0.003 


0.78 ±0.05 


550 


0.137 ±0.024 




0.157 ±0.005 




558 


0.117 ±0.010 






0.75 ±0.06 


584 


0.244 ± 0.060 


19.1 ±0.1 


0.108 ±0.003 


0.64 ±0.05 


600 


0.177 ±0.036 






0.43 ±0.20 


602 


0.052 ± 0.007 






1.76 ±0.10 


624 


0.077 ± 0.020 






1.30 ±0.05 


625 


0.197 ±0.058 




0.070 ± 0.003 




654 


0.043 ± 0.011 


20.5 ±0.5 


0.280 ±0.050 


1.46 ±0.10 


662 


0.193 ±0.028 






1.32 ±0.33 


674 


0.206 ± 0.033 






0.81 ±0.10 


704 


0.076 ± 0.010 


15.7 ±0.1 


0.305 ± 0.003 


1.45 ±0.10 


737 


0.136 ±0.043 






0.84 ±0.05 


787 


0.120 ±0.022 




0.087 ± 0.003 




796 


0.205 ± 0.041 




0.124 ±0.011 


0.98 ±0.02 


849 


0.115 ±0.016 






0.95 ±0.05 


857 


0.225 ± 0.026 






0.75 ±0.16 


863 


0.112 ±0.016 


18.1 ±0.2 


0.052 ±0.005 


0.40 ±0.10 


887 


0.230 ± 0.018 




0.101 ±0.050 


0.76 ±0.05 


925 


0.253 ± 0.053 


19.6 ±0.2 


0.065 ±0.005 




1036 


0.212 ±0.026 


20.6 ±0.2 


0.112 ±0.003 


0.84 ±0.02 


1052 


0.273 ± 0.074 






0.67 ±0.05 


1058 


0.242 ± 0.024 






0.69 ±0.10 
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Table 1:: (continued) 



Asteroid 


Pv 


a 


h 


1 min 


1105 


0.102 ±0.017 






1.20 ±0.16 


1355 


0.466 ± 0.082 


18.2 ±0.1 


0.083 ± 0.020 




1627 


0.153 ±0.046 




0.131 ±0.003 




1672 


0.094 ±0.016 




0.131 ±0.003 




1685 


0.292 ±0.127 




0.099 ± 0.003 




2131 


0.198 ±0.034 






0.86 ±0.15 


2577 


0.377 ± 0.062 


20.9 ±0.1 


0.124 ±0.044 




3169 


0.423 ± 0.067 


19.6 ±0.1 


0.276 ±0.018 




6249 


0.878 ±0.140 


22.4 ±0.1 


0.164 ±0.035 




6911 


0.454 ± 0.083 






0.83 ±0.16 



3. Revised Polarimetric- Albedo Relationship 



In Figures [2][4] we compare the measured WISE albedo to the slope of the polarization be- 
yond the inversion angle, the depth of the negative branch of polarization, and the inversion angle, 
respectively, for all objects with recorded values for these parameters. We distinguish objects in 
the Hungaria region and in the NEO population as red squares and cyan triangles, respectively. 
While the NEOs appear consistent with the MBAs, the Hungaria objects deviate from the general 
trend significantly. As discussed in lMasiero et al\ (|201ll ) the albedos for these objects are suspect: 
large deviations in the magnitude-phase slope parameter from the assumed G = 0.15 used for most 



asteroids can result in incorrect H values, and thus poorly constrained albedos ([Harris et aZ.l ll989) . 
Alternatively, large-amplitude long-period light curves may also corrupt the H values calculated 
from optical photometry. We are currently working on a program to better constrain the photo- 
metric parameters and albedos of these objects, but for the following discussion we will disregard 
the Hungaria asteroids. 

We see no overall trend between albedo and inversion angle in our data. The object with the 
anomalously high inversion angle is (234) Barbar a, the principal me mber of the "Barbarian" group 
of objects with strange polarimetric properties (jCellino et aZ.I I20Q6I ) . The object with an inversion 
angle well below the general trend is (704) Interamnia: some F-cla ss objects like Intera mnia have 
previously been shown to display unusually small inversion angles (jBelskaya et aZ.ll2Q05l ). 



We see the expected general trends when looking at slope and P m in, with low albedo objects 
showing steeper slopes and deeper troughs. We note, however, that the lowest albedo objects, those 
with pv < 0.04, have almost no representation in the polarimetric sample despite representing over 
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Fig. 2. — Albedo vs. slope of the polarization phase curve beyond the inversion angle. Objects 
located in the Hungaria region are noted as red squares, objects in the NEO population as cyan 
triangles. The green dashed line shows the best fit found for Equation [1] with our data. 

10% of the total population of MBAs observed by WISE, even before correcting the population 
for objects without optical followup (and thus without measured albedos). We therefore cannot 
comment on the reliability of the polarization-phase relation at the lowest albedos. A campaign 
of polarimetric observations of low albedo asteroids is critical to test these relations at their low 
albedo extreme. 

We find that the optimal description of the relation between albedo and polarimetric parame- 
ters is a linear fit in the three dimensional space of logpy-log /i-log P m in- We use only those objects 
where both polarimetric parameters are measured to an accuracy of 20% or better, leaving us with 
41 objects in our high-confidence sample. We perform orthogonal distance regression on the three 
dimensional data, using the associated errors on each measurement to determine the best fitting 
linear parameters as well as each parameter's error. We then reduce the best-fit parameters back to 
two-dimensional projections, which result in the following constant parameters for the relationships 
in Equations [1] and [2j 

Ci = -1.207 ±0.067 
C 2 = -1.892 ±0.141 
C 3 = -1.579 ±0.084 
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Fig. 3. — The same as Figure [2j but for albedo vs. depth of the minimum polarization branch. 
The green dashed line shows the best fit found for Equation [2] with our data. 



-0.880 ±0.106 



These projected fits are shown as green dashed lines in Figures [2] and E3 With the e xcept ion of 



C3, these parameters are all within one-sigma of the values found by iCellino et al\ ()1999l ). and 
all are within 1.5 sigma. As the WISE al bedos for the largest asteroids have been shown to be 
generally consistent with the IRAS values (jMainzer et adl201ld ) and all of the objects used here 
are in the size range sampled by IRAS, this agreement was not unexpected. Of the objects in our 
high-confidence polarimetric sample only 6 were not observed by IRAS, however the WISE albedos 
are all derived from a minimum of 5 observations (and an average of > 10) spread over time and 
thus are less sensitive to rotation effects. As the WISE data cover MBAs down to a few kilometers 
and NEOs to much smaller sizes, future polarimetric surveys focusing on smaller asteroids will 
allow this relationship to be tested over a more extensive size range. 

We can also project our fit of the polarimetric properties onto an axis of maximal variation. 
We define p* as the quantity of maximum polarimetric variation and find a best-fitting transform 
of: 

p* = (0.79 ± 0.02) log h + (0.61 ± 0.03) log P min 
Using our data we find a best-fitting relation between p* and albedo of: 



logp v = (-1.58 ± 0.09) + (-1.04 ± 0.04)p* 



-12- 




Fig. 4. — The same as Figure but for albedo vs. inversion angle. 



We show p* compared to albedo for all of the high-confidence objects with both measured h and 
P m in, along with this fit, in Figure We label (64) Angelina as " E", ( 2) Pallas as "B", and 



(132) Aethra as "M" following their Tholen taxonomic classifications (|Neesell2010l ). The clusters 
of objects with Tholen S and C taxonomic classifications are labeled as such, and include other 
objects within those taxonomic complexes (e.g. F- and D-types are included in the C-complex) . 
We indicate (71) Niobe with "n" in this plot; though it has a Tholen class of S it is distinct enough 
from the general cluster to warrant mention. Additional polarimetric and spectroscopic followup 
of this object will help determine why its properties differ from the general S complex. We focus 
on Tholen classifications he re as this taxonomic system shows the greatest distinction in albedo 
between the different types ([Mainzer et aZ.ll2Qlldl ). 



In addition to the lack of the lowest albedo objects in the polarimetric sample, we observe 
an over-representation of high albedo objects compared to the distribution for all similarly sized 
MBAs. Figure [6] shows the distribution of albedos of all objects in the sample with high-confidence 
polarimetric properties used to derive the linear three dimensional fit (the smallest of which is 
D ~ 35 km), as well as all MBAs larger than 30 km in diameter that were observed by WISE. The 
difference in these two distributions can be traced to the optical selection bias in the acquisition and 
measurement of the polarimetric properties of asteroids: very high signal-to-noise levels are needed 
to reach the polarimetric sensitivities that allow for accurate measurement of P m in and h. Thus, 
even though at a given phase angle low albedo objects will show larger degrees of polarization, the 
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Fig. 5. — The maximal variation of the polarimetric properties p* as defined in the text compared 
to the measured albedo. The best-fitting relationship projected from the three-dimensional fit is 
shown as the dashed line. Text labels denote specific asteroids or groupings as discussed in the 
text. 



reduction in photons received from these sources make these measurements less precise. Focusing 
future polarimetric surveys on low albedo asteroids will help make this sample more representative 
of the true distribution of MBAs. 



Conclusions 



Using the newly available albedos from the WISE space telescope thermal infrared survey data, 
we have fit the relationships between albedo, the slope of the polarimetric-phase curve beyond the 
inversion angle, and the maximum depth of the negative polarization trough. We restrict ourselves 
to objects with well characterized polarimetric properties (i.e. relative errors < 20%). Due to 
the selection of the polarimetrically observed objects this results in our sample consisting of only 
objects larger than D > 30 km, with nearly three-quarters having D > 50 km. 

We find that the function that best describes the albedo and polarimetry is a three dimensional 
linear fit in log py-log /i-log P m in space. Orthogonal distance regression allows us to find the best 
fitting parameters while accounting for measurement error on all parameters. When the best fit 
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Fig. 6. — Normalized albedo distribution of the high-confidence polarimetric sample (solid) and all 
MBAs larger than 30 km that were detected during the fully cryogenic portion of the NEOWISE 
survey (dotted). 



line is projecte d to two dime n sions , we find the resultant fit parameters are all within 1.5 sigma of 
those found bv lCellino et al\ (|1999l ). We also define a new polarimetric quantity p* that describes 
the maximum variation in polarimetric properties when compared with albedo. 

We observe distinct separation of some taxonomic classes in p* space. In particular E-type, 
B-type, and some M-type asteroids are far removed from the clumps that trace the more generic 
S- and C-complex objects. Asteroid (71) Niobe also holds a distinct location in this space despite 
its S-type classification under the Tholen system, and warrants further study. We note that the 
principal component (PC) analysis of Niobe from the Eight Color Asteroid Survey indicates that 
it is on the edge of the S-complex (jTholenl 11984 ) and it has a PC4 in the bottom two percent of 
all asteroids in th at survey (one of only 3 S-type or probable S-type objects with PC4 that low 
Zellner et Apod ). 



Finally, despite the prevalence of low albedo asteroids seen throughout the Main Belt (IMasiero et al 



201 ll ) we find that they are under-represented in the polarimetric sample. Notably, roughly 10% 
of MBAs have albedos py < 0.04, but there are no objects in our polarimetric sample with albedos 
this low. We recommend that future surveys focus on measuring polarization-phase curves for low 
albedo asteroids to properly sample this population. 
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